Impact of anisotropy on vortex clusters and their dynamics 
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We investigate the effects of anisotropy on the stability and dynamics of vortex cluster states 
which arise in Bose-Einstein condensates. Sufficiently strong anisotropics are shown to stabilize 
states with arbitrary numbers of vortices that are highly unstable in the isotropic limit. Conversely, 
anisotropy can be used to destabilize states which are stable in the isotropic limit. Near the linear 
limit, we identify the bifurcations of vortex states including their emergence from linear eigenstates, 
while in the strongly nonlinear limit, a particle-like description of the dynamics of the vortices in the 
anisotropic trap is developed. Both are in very good agreement with numerical results. Collective 
modes of stabilized many vortex cluster states are demonstrated. 



Introduction. Vortices constitute one of the nonlinear 
wave structures that have received considerable attention 
not only within Bose-Einstein condensates (BECs) in di- 
lute atomic gases [15 , but also in nonlinear optics [MS], 
among other fields [10]. Nevertheless, admittedly, BECs 
constitute one of the pristine settings where structural 
and dynamical properties of single- and multi- vortex or 
multi-charged-vortex states can be investigated experi- 
mentally and compared to theoretical predictions. It is 
for that reason that multiple techniques were developed 
to produce such vortices through phase-imprinting [TT] , 
stirring [12] or nonlinear interference [13] and few vortex 
states [14], vortex lattices [15 and higher charge states 
[T6] were generated experimentally. 

Although single- and multi-charge vortices were in- 
tensely studied, one of the themes that has received some- 
what less attention is that of clusters of few vortices. 
Central questions here include the possible geometries 
and structures of vortex clusters as well as their stabil- 
ity. Furthermore, the dynamics of cluster configurations 
is of immediate interest: due to the intricate interaction 
among the vortices, we indeed expect a rich "vibrational" 
dynamics of vortex clusters. In the literature mostly the 
vortex dipole has been considered [TTHSTj and it has been 
argued to be a robust dynamical configuration, while 
other states such as tripoles, quadrupoles, etc. are gener- 
ally considered to be fragile p~9j [21] . There are interest- 
ing connections between these states and dark solitonic 
stripes [20 - 22 , since some of the multi- vortex states arise 
from the instabilities of such dark-soliton stripes. The 
current interest in vortex clusters is also stimulated by 
several very recent experiments on the preparation of vor- 
tex dipoles [23, 24] and three- vortex states [25 as well as 
their dynamics. 

Our aim in this work is to present a detailed discus- 
sion of the existence, stability and dynamical properties 
of aligned vortex states by encompassing anisotropy i.e., 
we include the isotropic case of the above works as a 
particular limit and extend it to arbitrary anisotropics. 
In the case of co-rotating vortices (in BECs with rota- 



tion) , works in a similar spirit are those of [26 , although 
in the latter case the multi-vortex state may indeed be 
the ground state of the system. Here we focus on the 
anisotropic variants of the cases relevant to the above 
recent experiments. We employ two complementary the- 
oretical methods: chiefly and for large atom numbers, 
we extend the particle picture for vortices, developed for 
isotropic settings in [21], to the anisotropic regime. The 
resulting set of coupled ordinary differential equations 
(ODEs) is then used to gain insight into the stability 
of stationary multi- vortex configurations. On the other 
hand, we report on near-linear techniques valid in the 
limit of small atom numbers which allow us to make a 
very general prediction about the stabilization of a class 
of vortex states. These complementary theoretical ap- 
proaches are corroborated by the explicit numerical solu- 
tions of the corresponding 2D Gross-Pitaevskii equation 
(GPE). We identify the vortex clusters via fixed point 
methods for a wide range of atom numbers and strengths 
of the anisotropy. Subsequently a Bogoliubov-de Gennes 
(BdG) analysis of these states is performed, thereby 
yielding information on their linear stability. To exam- 
ine their robustness beyond linear stability theory, we 
disturb vortex states with white noise and employ real- 
time propagation. Our key conclusion is that anisotropy 
can controllably stabilize a whole class of unstable states 
of the isotropic limit, but it can also controllably desta- 
bilize stable states of that limit; these predictions should 
be directly testable in the experimental settings of refs. 
[23H25]. 

Setup. The general framework is set by the dimension- 
less GPE: 
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il>(x,y,t). 

In the BEC setting this implies that length, time, energy 
and density |^| 2 are measured in units of a Z: cj" 1 , hu z 
and (^l^f^hiaaz)^ , respectively; a is the s-wave scatter- 
ing length. uo z and a z refer to the oscillator frequency 
and length in the z-direction. The anisotropic trap is 



V(x,y) = (oo 2 x 2 + u 2 y 2 )/2. Stationary vortex states 
will be sought in the form ifj(x,y,t) = exp(—ijjLt)u(x,y), 
where \i denotes the chemical potential. 

Results. The single vortex, dipole and tripole states 
will be the prototypical ones considered here. They are 
shown as insets in Figs, [I] [2] and [3J respectively, for 
an anisotropic harmonic potential (with uj x = 0.2, uo y = 
0.08) and at \i = 2.5. We refer to these as "aligned 
vortex states". While the single vortex and the vortex 
dipole are stable in the isotropic limit, all higher aligned 
vortex states (tripole, quadrupole, quintopole, etc. [21]) 
are subject to dynamical instabilities in that limit. 

The vortex particle picture from [2TJ [27] can be nat- 
urally extended to take anisotropic confinement into ac- 
count. In the asymptotic regime of high values of /i, one 
can use the work of [28 to examine the consequences 
of anisotropy on the vortex precession. This yields the 
following equations of motion for a single vortex: x = 
—ujyQy and y = ou 2 Qx, where Q = ln(A/i/cj e ff)/(2/i), 
A « 2^/2tt is a numerical constant [2F and cj e ff = 

(uj 2 + ujy)/2. This implies that the precession fre- 
quency in the anisotropic trap is given by oo pr = uj x uj y Q, 
which yields the isotropic limit of [HJ HZ] when uj x = u y . 
Our numerical results on the single vortex state are in 
excellent agreement with the predictions from the above 
particle picture: At high enough values of /i, for every 
2D aspect ratio a = uo y /uo x that we have considered, the 
stationary single vortex solution exists, is dynamically 
stable and the value of the single anomalous mode in its 
BdG spectrum (responsible for vortex precession [27]) is 
accurately described by uo w above, as shown in Fig. [I] 
Note that in our numerical calculations we kept uo x = 0.2 
fixed and varied uo y to scan different values of a. Notice 
also that, in addition to the vortex internal modes, in 
both this and in following BdG plots, the BdG analysis 
captures the modes of the background state on top of 
which the vortex "lives". These include e.g. the dipo- 
lar frequency at uo = 0.2 due to uo x = 0.2, the dipolar 
frequency of oscillation around the y- axis (which is 0.2 
in the isotropic case, but varies as uj v varies), as well as 
higher order such as quadrupolar modes. 

If we now consider multiple vortices, then the corre- 
sponding large density particle picture should include the 
combined effects of the anisotropic precession and vortex 
interaction dynamics. Thus, the equations of motion of 
a vortex-antivortex pair i.e., a vortex dipole, will read: 

x 1 = -u y Qy 1 + B—^— (2) 

yx = u x Qx! - B ( 3 ) 

X2 = u y QV2 ~ B 2 ^ 2 (4) 

V2 = -u x Qx 2 + B ^ 2 , (5) 

where p 2 = (x\ — x^) 2 + (yi — y2) 2 and B « 1.95. These 
equations suggest that the two effects acting on the vor- 
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FIG. 1. BdG spectrum of the single vortex as a function of 
the anisotropy parameter a at \i — 2.5. The prediction for 
the precession frequency obtained from the particle picture 
is shown as solid black line which is almost completely cov- 
ered by the corresponding numerical data. The insets show 
contour plots of the vortex state's square modulus and phase 
profiles at a — 0.4, \i — 2.5. 

tices can balance each other in which case an equilib- 
rium position is obtained, described here by x\^ = 

and yi : 2 = ±<y/ B /(4uj 2 Q). One can also predict the lin- 
earized dynamics around this position, finding two pairs 
of frequencies: = iV^pr and uj^^ = ±u pT y/l — a 2 , 
where u pr is the precession frequency of a single vortex in 
the anisotropic trap. These two "anomalous modes" of 
the two-vortex state are shown in Fig. [2j together with 
the numerically calculated BdG spectrum of the vortex 
dipole. 
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FIG. 2. BdG spectrum of the vortex dipole as a function of 
the anisotropy parameter a at \i — 2.5; predictions from the 
particle picture are shown as black lines. The inset shows the 
dipole at a — 0.4, \i — 2.5. 

Here and in Figs, [l] and [3] the coloring of the BdG 
modes is chosen as follows: in the real part of the spec- 
trum, the anomalous modes (of negative energy) corre- 



sponding to intrinsic vortex motions are shown in red 
(gray), while the other modes which are modes of the 
background condensate are denoted in blue (light gray). 
The predictions from the linearization of the particle pic- 
ture ODEs are shown in black (dark). Notice again the 
excellent agreement between the theoretical prediction 
for the frequencies of the internal modes of the vortex 
pair and the corresponding numerical results. For any 
a 7^ 1 the rotational symmetry is broken and the corre- 
sponding mode in the BdG spectrum (which is present 
for any stationary solution in isotropic settings) deviates 
from its formerly vanishing eigenfrequency. Remarkably, 
for a > 1 this mode becomes purely imaginary, signal- 
ing destabilization of the dipole (and further destabiliza- 
tion of any higher order) configuration. This is in agree- 
ment with the ODE linearization results and with physi- 
cal intuition suggesting that anisotropy aligning the BEC 
with the vortex cluster should enhance its stability while 
anisotropy in the opposite direction should destabilize it. 

This particle approach can be generalized to any of the 
aligned vortex states. The case of three vortices ("vortex 
tripole") located in equilibrium at #1,2,3 = and 2/2 = 0, 

while y\ — —yz — -\Jb/(4:UJ 2 Q) leads naturally to three 

anomalous mode pairs (the number of such modes as- 
sociated with vortex motions is equal to the number of 
vortices in the cluster) with frequencies: 

cji 5 2 = ±V2u pr ^/l - a 2 (6) 
^3,4,5,6 = ±cj pr ^4 - 5a 2 ± y/9 + 2a 2 + 25a 4 . (7) 

In the isotropic limit of a = 1, these yield an unstable 
mode at iV^pr, a neutral one at (due to rotational 
invariance) and a precession frequency of y/Eu) pr . Con- 
cerning stability, the most important conclusion here is 
that there is a critical anisotropy of a = ~ 0.408, 

below which the tripole becomes completely stabilised. 
I.e., not only can anisotropy (in the "wrong direction" 
i.e., a > 1) destabilize stable vortex clusters such as the 
dipole, but also in the "right direction" (i.e., a < 1), it 
can stabilize not only the tripole, but in fact any higher 
order aligned vortex state (e.g. we have found stabiliza- 
tion in an aligned quadrupole not shown here). What 
changes is the precise critical value of a below which this 
stabilization occurs. Examining the numerically found 
BdG spectrum of the tripole shown in Fig. |3j we can 
see that the particle picture still captures the overall be- 
haviour of the relevant modes. Most importantly, stabi- 
lization of the tripole for small values of a is correctly 
predicted by the particle picture. The critical value of 
a obtained from the BdG analysis is given by a « 0.5, 
slightly deviating from the predicted critical anisotropy 
from the ODE approach. 

Stabilization of the aligned vortex states for small a 
can also be understood from the (near-linear) limit of 
the number of atoms N —> 0, in which the GPE re- 
duces to the linear Schrodinger equation with the chemi- 
cal potential \i playing the role of energy. In this regime, 
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FIG. 3. BdG spectrum of the vortex tripole as a function 
of the anisotropy parameter a at \i — 2.5; predictions from 
the particle picture shown as black lines. The inset shows the 
tripole state at a = 0.4, \i — 2.5. 



few-mode Galerkin-type methods have been developed 
in [2TJ [29] for the isotropic case which can be extended 
to arbitrary a. More specifically, for small particle num- 
bers the aligned vortex states and the soliton stripe states 
from which they bifurcate as the parameter \i is increased 
are well described by superpositions of harmonic oscil- 
lator eigenstates [21 . This picture is also applicable in 
anisotropic settings, and it is found that changing a (and 
thus shifting the energies of the linear eigenstates) quali- 
tatively affects the bifurcation diagram and, in turn, the 
stability of the states involved. 

In particular, the soliton stripe state with zero density 
along (and a phase jump of tt across) the ?/-axis can be 
traced back to the harmonic oscillator state ipio(x,y) = 
f ijji(x)ipo(y) in the linear limit, where ip n denotes the n- 
th excited state of the one-dimensional quantum har- 
monic oscillator. As \i is increased, the aligned vortex 
states subsequently bifurcate from this solitonic branch 
as has been shown for the isotropic case in Ref. [21]. 
For a = 1, except for the single vortex state all aligned 
vortex states bifurcate from the soliton stripe branch at 
finite values of the particle number. Since aligned n- 
vortex states emerge from the soliton stripe due to an 
admixture (to ^10) of ii^On oscillator components, and 
admixtures with high n are energetically suppressed, they 
occur at correspondingly high values of \i. In the isotropic 
case, the first vortex state bifurcating from the solitonic 
branch at finite N is the vortex dipole, then the tripole, 
quadrupole and so on. The dipole inherits the soliton's 
stability, while the soliton itself is destabilized due to this 
bifurcation. The tripole, quadrupole etc. then bifurcate 
from the already unstable solitonic branch, and thus are 
unstable themselves. 

In the following, let us apply the Galerkin-type ap- 
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proach (see details in [2TJ [29] ) to describe the bifurca- 
tions leading to the lowest aligned vortex states in the 
presence of anisotropy. The critical value of the parti- 
cle number and the chemical potential, respectively, at 
which the aligned n- vortex state is predicted to bifurcate 
from the soliton stripe branch are given by [21] 



Ma 



^0n ~ ^10 
: Wio - 



- B 



(8) 
(9) 



where A = j ifjf dxdy, B = j ipl Q ipQ n dxdy, and uiq = 
\uj x + \w y , cjon = \w x + (n+ \ )w y denote the energies of 
-010 and ^on, respectively. Note that in all cases relevant 
to our discussion A > B. 

Eq. (|8| implies that N cr crucially depends on the en- 
ergy difference of i/jio and ^o n ? which in turn is controlled 
by the value of the anisotropy parameter a. Choosing 
a > 1, the energy difference is increased, and the bifur- 
cations leading to aligned vortex states occur at higher 
values of N (or ja) than in the isotropic case. Physically 
speaking, the admixture of a i^o n component to ipio is 
energetically suppressed for large u y . As a special case 
of this, the degeneracy of ^oi and ipio underlying the 
single vortex state is lifted, and their superposition does 
no longer form a stationary state in the linear limit. In- 
stead, the single vortex bifurcates from the soliton stripe 
along the y-axis at a finite particle number. This is now 
the first bifurcation from the solitonic branch, leading to 
its destabilization. The dipole then bifurcates from the 
already destabilized soliton stripe and is thus unstable 
as well, in full agreement with the results from the parti- 
cle picture and our numerics. The single vortex branch, 
on the other hand, inherits the soliton's stability and is 
stable for any a > 1. The numerically determined bifur- 
cation diagram for a = 1.5 is shown in Fig. 4(a)[ together 
with the critical values of \i calculated using the Galerkin 
approach. 
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FIG. 4. Bifurcation diagrams for different values of the 
anisotropy parameter. Stable branches are shown as solid, 
unstable branches as dashed lines. Predictions for the bifur- 
cation points based on the Galerkin-type approach are given 
as dash-dotted vertical lines. 



Conversely, for values of a < 1 the energy difference in 
the numerator of Eq. ^ gets smaller, and the bifurcation 
points are shifted towards smaller values of TV and /i. 
Considering the single vortex state first, we find that it 
no longer bifurcates from the ipio soliton stripe branch in 
this anisotropic regime. This can be expected, as for a < 
1 the energy of ^oi is lower than that of ^io, and so the 
single vortex bifurcates from the energetically favourable 
-001 soliton stripe branch, characterized by zero density 
along the x-axis. As a consequence, the single vortex is no 
longer included in th e bifurcation diagram for a = 0.7, 
shown in Fig. |4(b)| Comparing to the a = 1.5 case 
clearly illustrates how the bifurcations of the dipole and 
tripole branches are shifted to smaller particle numbers. 

Lowering a even further continues this trend, and even- 
tually at a = 1/2 the vortex dipole arises in the linear 
limit, due to the degeneracy of ipio and ^02 • For values 
a < 1/2, the dipole bifurcates from the now energetically 
favorable two soliton branch starting from ^02 • Concern- 
ing stability, this means that the bifurcation of the dipole 
can no longer destabilize the soliton stripe branch for 
a < 1/2, which in turn implies that the tripole branch 
bifurcates from the stable soliton stripe branch and in- 
herits its stability in this regime. This, again, agrees well 
with the numerically found BdG spectrum in Fig. [3] 

The picture keeps repeating for decreasing values of 
a: Generally, for a = l/(n — 1), the aligned (n — 1) 
vortex state emerges in the linear limit, and the aligned 
n-vortex state takes its role in the bifurcation diagram 
and becomes stable, as it is the first state to bifurcate 
from the soliton stripe branch. This shows that the n- 
vortex state is stable in the regime 1/n < a < l/(n — 1). 
Additionally, for a < 1/n, the aligned n-vortex state 
is the first state bifurcating from the solitonic branch 
starting at ^o n , inheriting that state's stability. 

Thus, in total, we conclude that the state with n 
vortices aligned along the y-axis is stable for any a < 
l/(n — 1). In particular, aligned vortex states with an ar- 
bitrary number of vortices can be stabilized using strong 
enough transverse confinement. Note that this very gen- 
eral result agrees well with our numerical findings for the 
single vortex, vortex dipole and vortex tripole cases dis- 
cussed above. 

We corroborate the above general conclusions about 
the role of anisotropy and its different regimes with nu- 
merical simulations examining the robustness of the ob- 
tained aligned vortex states. To this end, a white noise 
signal is added to the state under consideration, and it 
is then propagated in time. Our observations are in ex- 
cellent agreement with the stability analysis: in regimes 
of a where the vortex states are found to be linearly sta- 
ble, they robustly persist. On the other hand, for a's 
where we find imaginary modes in the BdG spectrum, 
the perturbation leads far from the equilibrium. Specif- 
ically, four examples are considered in Fig. [5] The top 
panel confirms the persistence of the vortex dipole in an 
isotropic setting. The second row illustrates the desta- 
bilization of the same structure under a = 1.5 (with the 
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FIG. 5. Time propagation of aligned vortex states seeded 
with white noise for different values of a: (a) a = 1 dipole, 
(b) a = 1.5 dipole, (c) a = 1 tripole, (d) a = 0.4 tripole. 
All plots show |V>| 2 , elapsed time is given in the upper left 
corners, horizontal direction is x, vertical y. 



BEC elongating perpendicularly to the dipole axis) and 
its transformation into a dynamic state with one vor- 
tex weakly precessing close to the center, while the other 
counter-precesses far from it. The tripole dynamics is 
shown in the third and fourth rows: while the isotropic 
third row case suffers from instability, where again one 
of the vortices gets "expelled" from the BEC center, the 
anisotropic fourth row case of a = 0.4 restabilizes the 
tripole in a BEC elongated along the tripole axis. 

Finally, to illustrate that our picture of the aligned 
vortex states and their stabilization is valid even for very 
large numbers of vortices, we discuss a more exotic state 
with 17 vortices aligned along the y-axis. According to 
the general result reported above, this state should be 
stabilised for a < 1/16 = 0.0625, and indeed we find 
that for a = 0.06 it remains intact when disturbed with 
white noise. Remarkably, exciting the 17 vortex state 
with eigenvectors belonging to anomalous modes from its 
BdG spectrum leads to a behaviour strongly reminiscent 



of a (discrete variant of a) classical string along the y-axis 
with fixed endpoints, supporting a standing wave. Thus, 
the anisotropy transcending the dimensional instability 
barrier allows the observation of a fundamental vibration 
mode and the first few "overtones" as is shown in Fig. [6] 
Naturally, the period of such modes is directly associated 
with their BdG frequencies of our analysis. 
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FIG. 6. Time propagation of \i/j\ 2 of an aligned 17 vortex 
state when disturbed with eigenvectors belonging to different 
anomalous modes from the BdG spectrum (a = 0.06). One 
period of the fundamental vibration (a) at uj = 0.0609 and the 
first two overtones (b,c) at 0.0574 and 0.0539 are shown. Note 
the different length scales in the (horizontal) x- and (vertical) 
^-direction. 



Conclusions. In the present work, we demonstrated 
how anisotropy can be used to manipulate the stability 
and dynamics of vortex clusters. If used to elongate the 
system in a direction perpendicular to that of aligned 
vortex states, it generically destabilizes them. More im- 
portantly, as the system becomes increasingly elongated 
along the direction of the aligned vortices it can be con- 
trollably stabilized below a critical anisotropy and can 
exhibit even robust oscillatory dynamics reminiscent of 
classical strings. Critical values of the anisotropy for the 
stabilization of aligned n- vortex clusters are provided. 
We believe that both our near-linear bifurcation analysis 
and our particle approach including anisotropic preces- 
sions and vortex interactions can be of value to other 
fields where vortex dynamics is of interest, such as non- 
linear optics and condensed matter physics. 
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